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Abstract 

We discuss several issues regarding material homogeneity and strain compatibility for 
materially uniform thin elastic shells from the viewpoint of a 3-dimensional theory, with small 
thickness, as well as a 2-dimensional Cosserat theory. A relationship between inhomogeneity 
and incompatibility measures under the two descriptions is developed. More specifically, 
we obtain explicit forms of intrinsic dislocation density tensors characterising inhomogeneity 
of a dislocated Cosserat shell. We also formulate a system of governing equations for the 
residual stress field emerging out of strain incompatibilities which in turn are related to 
inhomogeneities. The equations are simplified for several cases under the Kirchhoff-Love 
assumption. 

Keywords: Material homogeneity, strain compatibility, shell theory, continuous distribution 
of defects. 


1 Introduction 

In this paper we explore the notions of material homogeneity and strain compatibility in materially 
uniform thin elastic shells. In general terms, a materially uniform body is said to be homogeneous 
if there exists a globally differentiable map from any conhguration of the body to its undistorted 
state P[I9]; this usually amounts to the body being free of topological defects. If a body is 
materially homogeneous then the associated strain held is compatible; the converse, however, is 
not true. Whenever material inhomogeneity leads to incompatibility, it subsequently becomes a 
source for internal stresses in the body [16]. A fundamental problem in micromechanics is, for a 
given distribution of defects in a solid, to determine the resulting state of deformation and stress 
held. The motivation for the present work is to develop a framework where this problem can be 
addressed for a broad class of thin structures. 

Homogeneity in materially uniform bodies has been explored both in 3-dimensional (3-d) solids 
[HiniiTSliTHliTHI and 2-dimensional (2-d) structured solids like shells [6l[8l[T2l[2Tl|23ll^|2^ but, 
unlike the former, appropriate inhomogeneity measures (or defect densities) and their relation 
with strain incompatibility has not been sufficiently developed for the latter. On the other hand, 
although strain compatibility relations for non-linear shells have appeared in the literature for over 
hve decades [HlEliTllilHlEn] , inclusion of incompatibilities has been attempted only recently M. 
Toward these ends one can model an inhomogeneous shell either as a 3-d body with small thickness, 
thereby using the available infrastructure of 3-d inhomogeneity theories, or alternatively as a 2- 
d Cosserat surface. It is well established that the 2-d approach is amenable to analytical and 
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numerical computations as well as for physical applications. It will therefore be useful to formulate 
inhomogeneity measures and incompatibility relations in terms of intrinsic surface quantities. Our 
methodology is to use previously established results from the 3-d theory to derive appropriate 
relations both in terms of 3-d as well as 2-d variables. In doing so, we need to dehne a 3-d metric 
field in terms of intrinsic 2-d strain measures; this dehnition essentially captures the geometry of 
the shell theory that we wish to work with, see Remark 13.21 for further details. A related notion 
of embedded homogeneity in thin structures, keeping in mind this dual characterization, has been 
recently discussed in the context of beams [TOlfTT] . Here, for a given a 1-dimensional materially 
inhomogeneous beam, the nature of material homogeneity of the 2-d solid, that it is actually made 
up of, is investigated. 

In Section 2, we revisit the classical fundamental theorem of Riemannian geometry in a new 
light, such that it is applicable to the 3-d continuum theory of topological defects. Several remarks 
are provided to clarify the usefulness of the theorem for material homogeneity, strain compatibility, 
and evaluation of residual stresses. Most of these results are well known both in differential geom¬ 
etry and in continuum theory of defects [3l|T6l[T9]. In Section 3, we apply the results from Section 
2 to explore the issue of homogeneity and compatibility in thin elastic shells. In order to use the 
3-d results, we construct a 3-d metric held using intrinsic strain measures of a Cosserat shell. More 
specihcally, we derive explicit relationships between 3-d continuous dislocation distribution and its 
2-d analogue on a Cosserat surface as well as representations of the latter in terms of Cosserat 
kinematical variables. We therefore obtain a complete characterization of the dislocation density 
distribution on a Cosserat shell. We discuss strain compatibility for shells, again within both 3-d 
and 2-d frameworks, and formulate the complete set of incompatibility helds while emphasizing 
their role in determining the residual stress held in the shell. We also derive the relationship of the 
incompatibility helds with intrinsic dislocation densities associated with the shell. In Section 4, 
we restrict our attention to Kirchhoh-Love shells and derive governing equations for residual stress 
determination under further simplihcations. In particular, we show that under Kirchhoh-Love 
constraint a dislocated shell can support only in-surface dislocations. 

2 Homogeneity and compatibility in a 3-d elastic solid 

Let H be a simply-connected open set in whose closure B has a piecewise smooth boundary; 
moreover, let B be such that it can be covered with a single chart. Hence, B admits a global 
parametrization X : {6^, 6^, 6^ =: (^) G —?■ B. The notation (•)^j is a shorthand for the partial 
derivative Although one can put various geometric structures (e.g., a connection, a metric) 
on B, it naturally inherits the Euclidean structure of including the Euclidean inner product 

(denoted by •)• Let Gi := X,i(r,C), Gij := Gi ■ Gj, := [Gij]-\ and G* := G^^Gj. The 
Roman indices vary between 1 and 3. 

We have the following 

Theorem 2.1 Let LF- he sufficiently smooth real functions on B, satisfying 

:= = 0 on (3. (2.1) 

(i) Then, there exists a sufficiently smooth invertible matrix field (denoted by [Hij{6°'X)]) on 
B such that 

L% = (H-y (2,2) 
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(ii) Moreover, let [gij] he a positive definite symmetric matrix field on B satisfying 

9ij;k • 9ij,k ^ik9pj ^jk9pi 0 (2-3) 

Then, 

9ij=9i-9j. (2.4) 

where Qj := HijG\ 

(Hi) Furthermore, if LT = then there exists a sufficiently smooth dijfeomorphism x ■ = X^G^ 
on B such that 

9ij = Xp\iXq\jG^'^, (2.5) 

where Xp\i ■= Xp,i - ^iXj, with TF ■= . + g^.^^ _ g^.^^). Here, [g^^] is the inverse of [gifi. 

This well-known result is a variant of the classical fundamental theorem of Riemannian ge¬ 
ometry (Theorem 1.6-1 in [3]) where otherwise the symmetry of LF is assumed a priori. The 
importance of asymmetric LF in describing continuous distribution of defects (or material inho¬ 
mogeneities) in a 3-d continuous solid body is discussed below. Several additional remarks are 
made before proceeding to the proof, so as to elaborate the utility of this theorem in describing 
the geometry of defects. In the next section, we will use this theorem to discuss analogous issues 
in a 2-d structured solid. 

Remark 2.1 (Materially uniform elastic solid) The trivial manifold B is our prototype for the 
theory of a continuous material body. The elements in B are called material points. The material 
structure of the body is modelled through a constitutive response function which can be used to 
understand the geometric nature of defects in the body. In the present article, we will assume the 
body to be simple hyper-elastic solid, for which the constitutive response function is given by a 
positive definite mapping W : Sym'^ x S —)■ M. Here, Sym~^ denotes the set of real, symmetric, 
positive definite matrices. In addition, the body is assumed to be materially uniform in the sense 
that there exist another positive definite mapping W : Sym~^ M and a matrix held [Hij] over B 
with det[Hij] > 0, such that 


W{IA„],X) = W([H„(X)][A„]{H„(X)]) (2.6) 

is satished for all [Aifi G Sym~^ and all X G H [19]. The held [Hij] is known as the material 
uniformity held and is, in general, not unique. 

Remark 2.2 (The material space) Let the manifold B be equipped with an affine connection with 
coefficients LF and a metric held with components g^j. The nature of the connection LF^ known 
as the material connection, and the metric gij, the material metric, is informed by the underlying 
material structure of the body in the following way. The geometry of the material space, dehned 
as the triple {B; LF,gij), brings out the defective nature of the material body. The Riemann- 
Christohel curvature of the material connection, as dehned by (I2.ip i . is a measure of the disclination 
content of the body. A zero disclination density (equation (12. ip ) is tantamount to the existence 
of a distant parallelism in the material space. This translates into the existence of well-dehned 
vector helds := HijG^ {[Hij] is the material uniformity held) which are covariantly constant 
with respect to the material structure (see fl2.17p L such that the material connection is necessarily 
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given by (I2.2p : are known as the material uniformity or crystallographic bases. On the other 

hand, the material metric Qij is derived from the nsnal Enclidean metric of the embedding space 
as a pnll-back of the material nniformity held; it is always well-dehned for a solid body [T9] . 
The non-metricity associated with the material space, dehned by the covariant derivative of the 
material metric with respect to the material connection (see (12.3^ 1 ). represents the presence of 
metric anomalies, snch as point defects or thermal strains, in the material body. Under zero 
non-metricity (eqnation (12.3p ). the metric is necessarily related to the crystallographic bases as in 
fl2.4p . The material space is, in general, non-Riemannian becanse the torsion tensor of the material 
connection, which has the components 


^ - rp. (2.7) 

is not necessarily zero. The third order skew tensor Th (or, eqnivalently, its second order axial 
tensor, with components := provides a measnre for the density of dislocation-like 

anomalies within the material body. The body is called materially homogeneous if and only if the 
Riemann-Christoffel cnrvatnre tensor, the non-metricity tensor, and the torsion tensor associated 
with the material space vanish identically at all points. 


Remark 2.3 (The Riemannian space) The metric strnctnre itself gives rise to another affine con¬ 
nection on B, i.e. the Levi-Civita connection with coefficients 

rfj := + 9qi,j ~ 9ij,q)j (2-8) 

which is, by dehnition, torsion free. It has however a non-zero Riemann-Christoffel cnrvatnre 


T^i _■pz _ pz _|_ p/z pz _ p/z pz 

'^jkl * ^ jl,k ^ jk,l ' ^ jl^ hk ^ jk^ hi"} 


(2.9) 


which provides a measnre for the incompatibility in the elastic Lagrangian strain field Eij := 
\{9ij — Gij). Thns, there exists a Riemannian space {B; gij^ associated with the material body. 


Remark 2.4 (Strain compatibility) Assnme Rh^^i = 0 and gij.^k = 0; hence the body is possibly 
dislocated. The material connection can then be shown to be related to the Levi-Civita connection 
as [19] 

= rf, + cp (2.10) 


where 

eg := - 9 '">’(rg, 9 ., + rp g„) ( 2 . 11 ) 

are the components of the contorsion tensor of the material connection Lh. Moreover, the respec¬ 
tive cnrvatnres are related to each other: 


Riti = + q-p - qn, + cgcL - (2.12) 

where the snbscript | denotes covariant derivative with respect to the connection TL. With Rt^^ = 0, 
the last relation provides a non-linear PDE for the metric for a given dislocation density held. For 

^In a Riemannian space, metric gij determines all the geometric structures. 
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a vanishing dislocation density, K^jki = 0 which is the compatibility relation for 3-d nonlinear 
elasticity. The inverse is however not true. A solution of the PDE (I2.12p with = 0 and 
^]ki = 0 is Ch = {Qy^Qij,i where [Qij] is a sufficiently smooth orthogonal matrix held over B 
(several other non-trivial solutions are given in [21]). When this is the case, the elastic strain is 
compatible even when the material is no longer homogeneous, in the sense that T?- ^ 0. Such 
a state of the material is known as contorted aeolotropy [I9|. For an isotropic elastic body in a 
state of contorted aeolotropy, however, vanishing of the curvature of its Riemannian space 
implies material homogeneity; this is due to the non-uniqueness in the torsion tensor for bodies 
with continuous symmetry groups [19]. Consequently, should be taken as a genuine measure 
of material inhomogeneity for an isotropic solid. 

Remark 2.5 (The residual stress held) The material space [B] L^-,gij) is the relaxed or stress- 
free state of the material body. Due to its general non-Riemannian nature, originating from the 
presence of material defects, it is often not physically realizable as a connected set in This in 
general leads to an incompatibility elastic strain held. For an elastic solid, in the absence of external 
forces and displacement boundary conditions, incompatibility of strain is the only source for non¬ 
trivial stress held [I6]. Equation (12.121) . with = 0, describes how dislocations as a source of 
material inhomogeneity yield strain incompatibility inside the material body and, together with 
the equation of motion and boundary conditions, form the governing equations for determining 
the stress held. 

We need the following lemma to prove Theorem 12.II (for a proof see Proposition 11.36 in [T7]b 


Lemma 2.1 Let V be a simply connected open set in M” and let [Aij{9^,z^)] be a sufficiently 
smooth [n x n) matrix field on V x M”, for (6*^,z^) eV x M"". If 


dAjj ^ ^ dAjj _ dAik , ^ dAik 
dQk Q^p ~ QQj P3 Q^p 


(2.13) 


onVx M", then given any (6*q, ..., 6*g) G V and {zl,..., Zq) G M", there exist unique smooth maps 
/i : 1/ —)■ R, for i = 1,... ,n, such that 


d f ■ 

^^=AffiO\h) onV, (2.14a) 

Mel,...,e-)=zl (2.14b) 


Proof of Theorem 2.1: (i) Let [H^] be an invertible matrix and let us consider the following 
partial diherential equation 

HijfiX) = L%{X)Hip{X), VX G B, (2.15a) 

HifiX,) = % (2.15b) 

for some generic Xq G B. For each integer / = 1, 2, 3, let fj := Hij and /° := Hiy Then, the above 
system of PDFs yield 

ffifiX) = L%{X)ffi{X), VX G B, 
fjiXo) = fl 
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(2.16a) 

(2.16b) 












The integrability condition of the above system is given by (I2.1jl . which follows from Lemma [2.11 
where Aij{9^,fm) := L'^{6^)fm- Equation (I2.15j) can also be written in the form (I2.2p . 

(a) Clearly, the vector helds ;= HijG^ satisfy the following problem by dehnition: 


:= - L%{X)g^{X) = 0, X e B, (2.17a) 

g^{Xo)=g^j, (2.17b) 

where g^ := i7°G*(Xo). Hence, the matrix held [g^ ■ gj] satishes 

(Qi ■ 9j),k = L^k i9p ■ 9j) + L^jk i9p ' 9i) in H, (2.18a) 

{9i-9j){Xo) = Qij, (2.18b) 


where := {Xo). The partial differential equation for gij given in the hypothesis (12.dh . 

along with the condition gij[Xo) = g^j, is identical to the problem fl2.18l) . The solution to the 
problem is however unique, see e.g. Theorem 1.6-1 in [3]. Hence, gtj = g^ ■ g^ in B. 

{in) When Lh = LL, the solution Hij to Equations (I2.15p satishes 

^li,j ^lj,i ^li\j '^61* ii' 

For a simply connected domain H, such that the curvature associated with Th is zero, Poincare’s 
lemma implies that, for each integer I (with values 1,2, and 3), there exists a sufficiently smooth 
Xi : H —)■ R such that 

Hu = xi\i in B, 

i.e., we have the existence of a smooth enough diheomorphism x fhat satishes (12.5p . □ 

3 Homogeneity and compatibility of thin elastic shells 

A shell is a 3-d solid body whose one of the dimensions is much smaller than the other two. 
Alternatively a shell can be described as a 2-d structured (Cosserat) solid. Our aim is to revisit 
Theorem 12.11 and the associated remarks, for a shell like body under certain assumptions on the 
nature of the deformation and the mechanical behaviour. We expect the resulting insights to be 
of value in studying distribution of defects, and the resulting stress held, in thin structures. Let 
us adapt the embedded curvilinear coordinates {9^,9'^X) ^ such that the coordinates {9^,9"^) 

he along the orientable mid-surface oo and ( along the normal direction to u. Then, B can be 
parametrized as 

x(r,c) = i?(n+cA^(n, (3.1) 

where R{9°‘) is the parametrization of u, N{9°') is the unit normal held on u, and ( G [—h, h], 
where 2h is the thickness of the shell assumed to be constant. The Greek indices take a value of 
either 1 or 2. Let Aa '■= R,a- The hrst and second fundamental form of oo are := Aq. ■ Ap 
and Baji := • Aq, respectively. Let [A“^] := and A" := A“^A^. With respect to the 

notation in the previous section, we have 

(i) GXe^X) = {€ - and ^ 3(^,0 = N{9X, where := 
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(ii) Gap = Aap - 2CBap + CCap and Gis = 5 * 3 , where Cap := BlB^p, 

(iii) + 2CA'^'^A^''B^^ + SCA^f^A^^C^^ + 0{C) and G*^ = 5 ^ 3 , and 

(iv) G“ = A“ + C G^ + 0{C) and = N, where G^ := A“^G*-^. 

Here, following Landau’s notation, for / ; M —)■ M^, we write /(s) = 0{s) if and only if there exist 
positive constants M and 5 such that ||/(s)||Rfc < M\s\ for all |s| < 5, and /(s) = o(s) if and only 
if lim.^o = 0 . 

3.1 Material homogeneity 

A motivation for our study is to derive intrinsic defect density (or inhomogeneity) measures for 
a 2-d structured body and obtain the required PDEs to solve for stresses and deformation of the 
inhomogeneous shell. Towards this end we start by considering the following set of strain-like 
smooth helds on uj\ Eap = Epa such that [Aap] + 2[Eap\ is positive dehnite, Aq, —1), Kap, 
and Aq. The choice of these strain-measures, whose physical nature is remarked below, dictates 
the kinematical nature of the thin structure. For instance, the straining of a membrane can be 
described in terms of Eap only. 


Theorem 3.1 (i) Let be sufficiently smooth real functions defined onB = \X(9°‘, () = i?(6*“)-|- 
CAf(r) C C e [-h,h], satisfying 



i?},, := L],^, - L],^, + L%E^, - L%L\, = 0 on B. 

(3.2) 

Let [gij] he a 

symmetric matrix field on B, defined by 



QoifS • ^a/3 "t“ C PotjS C Qa^-} QcuS • C fi^33 • ^5 

(3.3) 

where 

^a/3 • -^ajS 

(3.4a) 



(3.4b) 


Qa0 • ^ 

(3.4c) 


Ua := — B^a) + Aq,(A + 1 ), and 

(3.4d) 


l/:=a“^A„A^-f(A + l)2, 

(3.4e) 

with [a"^] := 

[(^ap\~^ , such tha^ 



9ij-,k ■ 9ij,k ^ik9pj ^jk9pi 0 

(3.5) 


^Note that, since uj is bounded and [gij\ is continuous in 0“ and C, [pif will be positive definite on := ux (—h, h) 
for sufficiently small h. Our result is valid for this sufficiently small h and we start with an B such that h conforms 
with this small value. For a technical discussion on the issue of smallness of h and positive definiteness of \gij\, 
please refer to the proof of Theorem 2.8-1 in [3]. 
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Then there exist sufficiently smooth vector fields aa{0'^) with cii( 6 *“) x 02 ( 6 *") 7 ^ 0, Da{0'^), and 
d{6°') with d[9°') ■ ai( 6 *“) x 02 ( 6 *“) 7 ^ 0 , such that (here should be interpreted as introduced in 
the previous section) 


C) = aa( 6 '") + and (3.6a) 

gfiefiC) = die^). (3.6b) 

In particular, 

ao( 6 >“) ■ a; 3 ( 6 >“) = a„/ 3 ( 6 >“), (3.7a) 

Da{9°‘) = (Ao-q, — + A^n, and (3.7b) 

d{9°') = Ao-d®" + (A + l)n, (3.7c) 


where a‘^(9°‘) := and n( 0 ") := ^ I, ■ 

^ ^ ^ ^ ^ |ai(0“) X a2(^“)| 

Furthermore, if the torsion tensor evaluated at the mid-surface vanishes, i.e. Tj^ (6*",0) = 0, 
then there exists a sufficiently smooth diffeomorphic image Cj of 00 , parametrized by r(6*“), such that 
«'a(6*") = r^a{9°‘), and £>q,(6 *") = d^a{9°‘). In particular, there exists an open sufficiently smooth 
diffeomorphic image B of B, parametrized by x{9°'X) = + fdifiX, such that gfid'^X) = 

A,ir,c)- 


This result is central to our understanding of the defective nature of the shell. We would like to 
emphasize that the hypothesis fl3.3l) for the metric structure on the 3-d body manifold B of the shell 
is crucial to our proof as well as for interpretation of the results. Relation fl3.3p is a generalization 
of a similar 3-d metric used in the proof of the fundamental theorem (Theorem 2.8-1 in [3]) of 
embedded surfaces in M^. Before proving the theorem, we provide several remarks to bring out 
the importance of this result from the view point of defect mechanics. 


Remark 3.1 (Shell kinematics) Consider a sufficiently thin shell (i.e. h/R « 1, with R = min¬ 
imum principal radius of curvature of the shell mid-surface for a given deformation) made up of 
an hyper-elastic solid whose material structure is characterized by a material connection Lf- and a 
metric gij, such that the curvature associated with LA is zero (see fl3.2p ) and the metric is covari- 
antly constant with respect to the material connection (see (13.51) ). Therefore the shell material has 
no intrinsic disclination and metric anomalies. Furthermore if this shell is free of dislocations, i.e. 
if C) = 0 , it will have a coherent relaxed conhguration characterized by a diffeomorphism x 

on B. The notion of sufficient thinness of the shell is then manifested in the particular action of x 
which maps 

B = {X{9^X) = RXX + CNfiX} (3.8) 

onto the relaxed state 

B = xiB) = {x{9fi 0 = rfiX + CdfiX}. (3.9) 

Here, the stress relaxation process (i.e. the elastic deformation) respects the classical Green-Naghdi 
approximation for sufficiently thin shells [13]. The elastic deformation x '■ B B of a materially 
homogeneous sufficiently thin shell preserves the ‘hbrous’ structure of B, in the sense that straight 










transverse sections remain straight and transverse throughout the deformation. Clearly, there are 
two separate modes of deformation at play: the mid-surface deformation 

i^(r) ^ r(r) = i^(r) + w(r), (3.10) 

where u{6°‘) is a well-dehned displacement held on u, and the director deformation 

N{e^) ^ d{e^) = Q(r)Af(r), (3.11) 

where ;= cZ(0") C) A/’(6*“) is a second-order tensor held dehned on u. 

Thus, while as a full-hedged 3-d material body the elastic deformation of a materially homoge¬ 
neous shell is characterized by the diheomorphism X (0“, Q i-^ x{6°', (), due to its ‘hbrous’ nature— 
a mathematical artifact brought forth by the representations (I3.8jl and (13.91) — the same elastic 
deformation is equivalently characterized by the set consisting of a diheomorphism It{9°‘) i—>■ r{9°‘) 
of the base manifold uj and a linear isomorphism N{9°‘) i-4 d{9°‘) of the transverse hbres (the 
director helds) attached to the base manifold. The latter characterization is a manifestation of the 
Cosserat structure of the shell described by a vector bundle consisting of a 2-d base manifold u 
and an isomorphic copy of a frame of attached to each point of G ca, such that two vectors 
of the frame at R span Tru and the third vector is transverse to it. In the Cosserat picture, the 
elastic deformation of a materially homogeneous shell is a principal bundle isomorphism of this 
vector bundle which subsumes a diffeomorphism of the base manifold and a linear isomorphism of 
the frame held, as just described!! 

Remark 3.2 (Strain measures) The dual characterization of shell kinematics, as outlined above, 
motivates the particular form (13.31) of the metric structure on the 3-d shell manifold B, relating 
the two equivalent descriptions of the shell material space. The dehnition (13.3p is central to all 
the subsequent claims made in Theorem 13.11 We are considering a single director geometrically 
nonlinear shell theory, taking into account the transverse shear and normal deformation. The helds 
«.a/ 3 , Aq/ 3 , Aq Aq, and A constitute the 2-d metric structure on the shell’s Cosserat material space. 
‘Sufficient thinness’ of the shell is encoded in the dehnition (13.3p which describes how the 2-d strain 
helds can be used to construct a 3-d metric on B. The resulting 3-d metric is second-order in the 
transverse coordinate C; it is the unique generalization of the metric 


C) = aa/3 - 2C&a/3 + (3.12) 

dehned in the proof of Theorem 2.8-1 in [3] that considers a Kirchhoh-Love shell model (i.e., 
Aq, = Aq = a = 0), where hap := — Aq^ -f Bap. Due to the specihc form (13.3p of Qij, the material 
uniformity bases on the 3-d material space {B; L^j,gij) decompose into the form (13.6p . guaran¬ 
teeing the existence of the vector helds d and that constitute the material uniformity bases 
on the Cosserat material space of the shell. The relations (13.7p . which are analogous to Theorem 
\2.Vi ii). relate the metric structure {oq^, Aq^, Aq, Aq, A} to the material structure {a^^Da^d} in 

point worthwhile overemphasizing is that nnlike the inherently Cosserat-type materials (also known as polar 
media or materials with microstructures) e.g. liquid crystals, magnetic materials etc., a structural shell (as is 
presently the case) is not inherently a Cosserat material but merely a mathematical artifact originated from the 
representations (13.811 and (13.911 . 
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the Cosserat material space of the shell. The last two relations in (13.Th can be more explicitly 
written as 


Aa = d ■ aa — N ■ Aa, (3.13a) 

A = d-n- N ■ N, (3.13b) 

Aaj 3 = Dp ■ tta — N^p ■ Aa, and (3.13c) 

Aa = Da ■ n — N^a ' N- (3.13d) 


These expressions clearly bring out the strain-like nature of the helds appearing on the left hand 
side. 

Remark 3.3 (Material inhomogeneity measures) Let us consider the 3-d shell to be free of disclina- 
tions and metric defects. According to Theorem I3.1f ii). for the kind of shell presently considered 
(characterized by a metric of the form fl3.3p L the material inhomogeneity of the shell is completely 
determined by the restriction of the 3-d torsion tensor to the mid-surface, i.e. Tk{d°‘,0). This 
tensor held can be alternatively expressed in terms of the 2 -d inhomogeneity measures of the shell: 


0) = 0)T^„^(r) + 0)T,apiO-), (3.14a) 

0) = (iL-i)3-(r, 0)T^„^(r) + (^-')33(r, 0)T3„^(r), (3.14b) 

0) = 0)T^a3{0n + O)T 3 „ 3 ( 0 “), and (3.14c) 

0) = O)T^„ 3 ( 0 “) + O)T 3 „ 3 ( 0 “), where (3.14d) 

T;,a/ 3 ( 6 '“) := 0) = (3.15a) 

T 3 „^(r) := 4[„,^](r,0) = i/ 3 [a,/ 3 ], (3.15b) 

T^o 3 (^°') •= -^/i[o, 3 ](^°') 0) = Fii3,a — Ffia + and (3.15c) 

T3«3(^“) := i^3[a,3](0“, 0) = F33,„ - ^ 3 ^, (3.15d) 

are the 2-d Cosserat inhomogeneity measures, with Hap{9°‘) := Aa ■ ap, ■= ■ aa, 


Fap := Aa ■ Dp, Fsp := N ■ Dp, Fas := Aa ■ d and F 33 := N ■ d. In the above expressions a square 
bracket in the subscript indicates the skew part of the held with respect to the enclosed indices 
(on the other hand, round brackets are used to indicate the symmetric part). The vanishing 
of 2-d Cosserat inhomogeneity measures is equivalent to T[j{9°‘,0) = 0. We note the following 
interpretations: 

1. The component T^q,^( 0") measures in-surface dislocation density. Its vanishing implies the 
existence of the in-surface components ra{9°‘) := r-Aa of the surface diheomorphism x(^", 0 ). 

2. The component T 3 q,^( 6 '“) measures the out-of-surface dislocation density. Its vanishing im¬ 
plies the existence of the out-of-surface component r^{9°‘) := r ■ N of the surface diheomor¬ 
phism x( 0 ", 0 ). 

3. The component T^o 3 ( 0 “) measures the in-surface integrability of the director held d. Its 

vanishing implies {Da — d^a) ■ = 0 . 
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4. The component T 3 q 3 ( 0 “) measures the out-of-surface integrability of the director held d. Its 
vanishing implies {Da — d^a) ■ ■«■ = 0. 

In the Cosserat picture the material inhomogeneity which encodes the integrability 

of the surface material uniformity bases should be interpreted as a density of dislocations 

smeared over the base manifold ca. On the other hand, the material inhomogeneity TjQ, 3 ( 6 *") should 
be interpreted as an ‘apparent’ disclination density, smeared over the base manifold ca, for they 
encode the compatibility of the out-of-surface material uniformity bases d{6°‘) and Da{0°') in the 
sense that a certain derivative of the transverse director held dehned as Vod := d a — Da, 

vanishes if and only if TjQ, 3 ( 0 “) = 00 The diherential operator Vq, gives rise to a parallelism, 
hence an Ehresmann connection [7] , on the normal sub-bundle of the original vector bundle in the 
Cosserat picture. We can summarize the above result as 

Proposition 3.1 Under the assumptions made in the present remark, the shell is materially ho¬ 
mogeneous if and only ifTiaj{9^) vanish simultaneously. 

Remark 3.4 Our framework generalizes the earlier notion of material homogeneity, and the asso¬ 
ciated measures, in a thin shell as was proposed in [ZllH]. The notion of material homogeneity 
discussed therein is the existence of a globally hat diheomorphic conhguration u, carrying every¬ 
where a normal director attached to cu. As is clearly evident from the statements of Theorems 
12.11 and 13.11 and the ensuing discussion, this should not necessarily be the general case because 
the diheomorphic image of the current conhguration of the shell (or, in the Cosserat picture, the 
image of the current vector bundle under a principal bundle isomorphism), upon stress relaxation, 
may not be globally hat at all. 

Proof of Theorem 3.1: According to Theorem 12.11 hypothesis (13.2p implies that there exist 

vector helds •= HijG^, where [Hij{6°‘,Q] is an invertible matrix held on B, satisfying 

fl2.2p . Let gi{0'^,C) be analytic in (. Then there exist vector helds on cu, aa{6°‘) := g'„(6*",0), 
Da{d°‘) := ga{0°',0), d{9°‘) := g^{9°‘,0), and Ea{9°') := 913 ( 6 *“, 0) (superscript prime denotes a 
derivative with respect to C), such that 

pjr, C) = aa{9^) + CDa{9n + o{C) and 

^3(r,c) = d(n+c^^(n+o(c). 

Hypothesis (13.51) subsequently implies that gij{9°‘,Cf) = gi{9°‘,C) ' 9j{9°‘,C)- From ga/3{9°',C) = 
C) ■ 9g{9'", C), we obtain 

aay{n + CPay{9n + C"Qag{9n = a„(0“) • a^(r) 

+C(a„(r) ■ D^{9^) + a^(r) ■ E>„(r)^ 
+eDa{9^)-Df,{9-) + 0{C^), (3.16) 

^ We interpret the quantity Vq< 1(0“) as an ‘apparent’ disclination density because this quantity contributes 
nothing to the curvature of the material space of the actual 3-d shell (while it does contribute to the curvature 
Kjj^i of its Riemannian space, hence producing incompatibility in the strain field). Only in the Cosserat picture, 
which is a mathematical artifact in the present case, a part of the actual 3-d dislocation density ‘appears as’ a 
disclination density. 
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which implies that g^{6'^,() is linear in (: 


gjr,c) = a„r) + ci^ar)- 

From C) = C) ■ C), we obtain 

A„ + C{a"^A,(A^„-5^„)+A„(A + l)} = a„(r) • d(r) 

+c(£>„(r) ■ d(r) + a„(r) ■ Eie^)^ 
+0{e), (3.17) 

which implies, hrstly, that fi'o 3 ( 6 *“,C) is linear in (: 

53(r,c) = d(n+c^^(n- 

But, since 5-33(6'", C) = fl'3(6'", C) ' fl'3(6'", C), i-e. 

A„A^ + (A + 1)2 = d(r) ■ cZ(r) + 0(C), (3.18) 

we obtain E{9°') = 0 or that 9 - 3 ( 6 *“, C) is independent of C,'. 

93(r,c) = d(n- 

Finally, relations (13.71) can be obtained using fl3.16p . (13.171) . and (I3.18p . 

[a) We can write various components of the torsion tensor, in terms of [Hij] and its inverse 
as 

(3.19a) 

+ (i7-i)33^3K/3], (3.19b) 

+ {H-^r^H 3 [a, 3 ], and (3.19c) 

= (i7-i)3“.ff,[„,3] + (F-1)33^3[„,3]. (3.19d) 

IfTf^(r,0) = 0, i.e. 77,;/ 3 ]|^_o = 0 (the covariant curl will then be same as the ordinary curl), 
we have fl'/ 3 ,Q,(^", 0) = ^o,,^(6*“, 0), which implies that there exist three sufficiently smooth real 
functions ri{9°‘) such that Hia{9°‘,0) = rj_Q,( 6 *“) and 

^jr,o) = a„(n = ^,«(n, 


where r = raA°‘ + r^N. 

Moreover, if T^^(9°‘,0) = 0, i.e. .^i[a; 3 ] = 0, we have d^ai9°‘) = 9q,, 3(6*“,0) or equivalently 

d„(r) = r>„(r). 

Evidently, we have the existence of a sufficiently smooth diffeomorphism x( 0 “, C) = r{9°‘) +C< 7 ( 6 *“) 
on B such that C) = X,i(^", 0- □ 
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3.2 Strain compatibility 

For a materially homogeneous thin shell, the components of the Riemann-Christoffel curvature 
of its Riemannian space [B] gij) vanish giving rise to the strain compatibility conditions. In 
this subsection, we will show that the strain compatibility relations for the 2-d strain measures 
can be recovered from the vanishing of Kijki at uj alone, not necessarily on the full B. A more 
restricted version of the following theorem has appeared in the proof of Theorem 2.8-1 in [3]. 

Theorem 3.2 If the restriction of the Riemannian curvature of the shell to ( = 0 is zero, i.e. 
z/0) = 0, then there exists an open sufficiently smooth diffeomorphic image B of B, 
parametrized by x(6*", () = r(9°‘) -|- such that gi{9^, C) = 0- 

Remark 3.5 (2-d strain compatibility conditions) For a materially homogeneous thin elastic shell, 
the intrinsic strain measures are given as 



(3.20a) 

5 

<! 

(3.20b) 

A = d - 1 , 

(3.20c) 

■^a/3 -|- da\l3 

(3.20d) 

^a. dfj 

(3.20e) 


where d = daOR + dn and aap = Viewed as a set of partial differential equations for r(6*") 

and d{9°‘), which together represent a directed surface u, the integrability conditions of the system 
(13.201) give us the compatibility conditions for the 2-d strain measures Eay, Aa/s, Aq, Aq, and A of a 
sufficiently thin shell. We now discuss these conditions, following Epstein |6j, before moving on to 
the proof of the theorem. The question of compatibility can be posed as: given sufficiently smooth 
helds including a symmetric matrix [Eay], an invertible matrix [A^/j], two vectors {A„} and {Aq,}, 
and a scalar A{^ —1) on a parametrized surface u, with its hrst and second fundamental form as 
Aais{9°‘) and Bay{9°‘), respectively, what are the conditions to be satished by the five given helds 
for the existence of a sufhciently smooth parametrized surface A, along with a director held d{9°‘), 
having its hrst fundamental form given by Oa/s = -l- 2 £'q,^ and its second fundamental form 

bay suitably constructed out of the given helds, so that the equations fl3.20p are satished. Towards 
this end, let us hrst ensure the existence of surface A. We can read oh the formula for its second 
fundamental form bap from the dehnitions fl3.20bl) . fl3.20bl) . and fl3.20cl) as 


^ajS 


AajS BafS Aa\y 

A + 1 


(3.21) 


To ensure that the given helds indeed give rise, via (I3.2ip . to an admissible second fundamental 
form for a realizable surface, we must ensure that biay] = 0, or equivalently 

J := e°'^ba0 = 0 or A[ay] - A[a\y] = 0, (3.22) 

where = Cay is the 2-d permutation symbol. Equation (l3.22p o is the hrst strain compatibility 
condition. Keeping fl3.22p in mind, the conditions on Oay = Aay + 2Eay and h(^ay) so that they 
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indeed constitute the first and the second fundamental forms of a parametrized surface ca, upto 
isometries of are the Gauss and Codazzi-Mainardi relations: 


STafia ^(at7)^(/3T) ^(o/?) ^(cr) aud 

b{aa),l3 — b(a/3),a + ~ ~ 0 ) 

provided u is simply connected, where 


(3.23) 

(3.24) 




TOi^CT 


\= s. 


ao-T,^ 




— 




^affi ■ 2 (®Q:/x,/3 d" aud 

^af3 •” ® ^^afii- 

Equations fl3.23p and fl3.24p are the second and the third strain compatibility conditions. The 
three compatibility conditions derived so far ensure existence of a unique (modulo an isometry in 
M^) surface u. Finally, to ensure the existence of a director held d(0“) on Cj, consistent with fl3.20p . 
we require 

If, := - A,f, + ^ 0 _ (3 25 ) 

This is the fourth compatibility condition, obtained by eliminating bfaf) (recall that = 0) from 
equations fl3.20bp - fl3.20el) . Altogether, we have six independent strain compatibility relations for 
twelve independent components of strain measures. 


Remark 3.6 (Residual stress held) Given the strain energy density W{Eij) per unit volume of the 
shell material, the 2-d strain energy density per unit area on ca, denoted by Aap, Aq, Aq, A) 

(or equivalently as U{aa,Da,d)), can be calculated via the following integration [13] 

ilj{Ea/3, Aaf, Aa, A^, A) = j VgW (Eij) (!(. (3.26) 

Then, with zero body force, the 2-d equilibrium equations for various ‘stress’ measures of the shell 
are 


T“ = 0, (3.27a) 

M" - fc = 0, and (3.27b) 

tta X T“ -I- Da X M°‘ + d X k = 0, (3.27c) 

where jT“ = da^U, = Od^U, jk = dJJ, j := a := det{aap), A := det{Aa 0 ), and 

G := det{Gij). If in addition the shell is materially inhomogeneous then the compatibility relations 
are not satished and provide additional equations for determination of the stresses, provided the 
incompatibility is known. We introduce six incompatibility measures J, K, Lg-, and such that 


the strain compatibility equations take the form 

J = J(r), (3.28a) 

*S'i2i2 ~ ^11^22 + ^q2) = IK(^“), (3.28b) 

^(f7i)|2 - ^<72)11 = IL^(6'"), and (3.28c) 

la = la(0“). (3.28d) 


14 














The incompatibilities are related to curvature fields K 1212 , Kuas, and according to the 

following six relations: 

K(r)=i^i2i2(0“,O), (3.29) 

+ (A + l)Li = Ki213{9^, 0), (3.30) 

- a^^A^K + (A + 1)L2 = i^i223(0“, 0), and (3.31) 

(A + 1) I(p|o-) — A(plo-) — CL J “ 1/3 j- 

+(A + l)j(^o“^ + (^(a + 

+6up(A^o — Bj3(r) + e^o(Ac.p — £?up)l — 0). (3.32) 

The above construction is analogous to Kroner’s framework [16] of residual stress determination 
for a given incompatibility tensor field in the context of 3-d linear elasticity. The incompatibilities 
can also be obtained from defect densities. For instance, consider the case when only dislocation 
anomalies are present in the 3-d shell. Then incompatibilities can be written in terms of the torsion 
tensor as 

K(9“) = -gpi(9",0)|c''2|, - + Cic;, - ( 3 , 33 ) 

+ (A + 1)L, = -9p,(9“, 0)lC'3n - C',,, + C", - (3.34) 

- o‘'’AjK+ (A + 1)L2 = -jp,(r,0)|Cf3p - C '213 + and (3.35) 


(A + l)I(p|p) A(plp) a ^ 

+(A + l)j(^a“a + 

(-^/3f7 B “h 6y3(j(-^ a.p 

= -g„(e‘, 0 ) IC',!, - ca ,|3 + CjjCJ, - clci,] I 

Proof of Theorem 3.2: The coefficients of the Levi-Civita connection, defined by T^p 
9jp,i 9ij,p )) 


Tsss — 0 , Tssp — Up — ^sps — Tpsa — 


bspcr Fp 3 (j 


A[o-,p] T ^ A Qpcr^ ) 


rp(T3 ^(a,p) n^pgr T C ( U(cr p) Qpa ) ) aud 


c 


FpcriS SpcrS T I PpS,a T Pa 5 ,p Pap,S 1 „ I QpS,a T Qad,p Qap ,5 


e 


(3.36) 


1 

2 


idipj T 


(3.37a) 

(3.37b) 

(3.37c) 

(3.37d) 
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The curvature Kijki has six independent components such that Kijki = 0 if and only if K 1212 = 0, 
Ki 2 a 3 = 0, and = 0. After some manipulations, it can be shown that 

-^i212|^^q := (r22i,i ~r2ii,2 + r2ir2ii ~ r22riii) 

= >S'i 212 — \bllb 22 — ^( 12 )]- (3.38) 

Hence, A'i 212 |^^q = 0 implies that 

<S'l 212 = ^ 11^22 — ^( 12 ); (3.39) 

which is the single independent Gauss’ equation, cf. (13.231) . Moreover, 

Kl2a‘i\(^^Q '■= (r231,cr — r2o-l,3 + r^o-Tsii — r 23 ro-li) I 

= a“^A/3 ^S'q,o-21 + \b[a2)b{al) — ^((Tl)^(a2)] ^ + (^ + 1) [^(o-l)|2 ~ ^(fT2)|l] , (3.40) 
which can be used to calculate 

A"i213|^^q = ^>S' 2 i 2 i + [&(i 2 )^( 2 i) ~ ^(ii)^( 22 )] ^ + (A + 1) [ 6 (ii )|2 — &(i 2 )|i] and (3.41) 

A"i223|^^q = —n^^A^ ^S'i212 + [^(21)^(12) ~ ^(11)^(22)]^ + (A + 1) [6(21)|2 ~ &(22)|l] • (3.42) 

Substituting fl3.39l) in fl3.41l) and fl3.42p . along with Ki 2 uz\^_q = 0, yield 

^(ii )|2 ~ ^(i 2 )|i = 0 and &( 2 i )|2 ~ &( 22 )|i = 0, (3.43) 

which are two independent Codazzi-Mainardi equations, cf. (I3.24p . Finally, we consider 

A'p3o-3|^^0 (r33p,(T - + Fg^Fspj - rggFo-pi) 

= (^ + 1) -^(pk) — ^(p-^o-) ~ n + 2^;fl(p'^o-) ~ '^'^1 

~ ^0(r) + ~ ^ap)'^ ■ (3-44) 

With A"p 3 o- 3 |^_Q = 0, 03.441) is a set of three coupled hrst order homogeneous non-linear partial 
differential algebraic equations (PDAE) for three unknowns and J. Note that 

(a) Jq, = 0 and J = 0 are solutions to this PDAE, but there can be other non-zero solutions 
whose nature depends strongly on the coefficient functions. 

{b) If we assume J = 0, the PDAE reduce down to a hrst order homogeneous overdetermined 
system of linear PDEs in I a- The system has a zero solution and other non-zero solutions depending 
on the coefficient functions. We disregard the unphysical non-zero solutions because they become 
unbounded under generic perturbations of the initial condition for generic coefficient functions [ 2 ]. 

(c) If we assume = 0, the PDAE reduce down to three quadratic algebraic equations in J 
which can be easily shown to have the unique solution J = 0. 

These three facts imply that /„ must be proportional to J, i.e. there exist functions £q,(J; 0") 
such that la = Ca{J',0'^)J. Using this in the PDAE reduces it to a set of overdetermined hrst 
order non-linear PDEs in J which clearly has a zero solution along with other unphysical non-zero 
solutions. The zero solution implies /„ = 0. □ 
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4 Kirchhoff-Love shell with continuous distribution of dis¬ 
locations 

We now restrict ourselves to the case when the inhomogeneous shell has only dislocation anomalies 
and the stress relaxation process respects the Kirchhoff-Love constraint d = n, i.e. Aq, = A = 
A„ = 0. As a result, the Cosserat material uniformity basis determines the complete Cosserat 
material uniformity bases. The compatibility relations /„ = 0 are trivially satished; hence, 
cannot take a non-zero value. In addition, the components (or equivalently Tjgo.) of the torsion 
tensor are identically zero as the director helds are compatible in the sense that Vo-n = 0. Under 
the present simplihcation, (13.361) reduces to 

Qj ^Jf Cap&jSa T -|- C/ 3 ai,-^ap 0. (4.1) 

This relation implies that either (j CapC^^ + eap{A.p^ — + ep„{Aap — Bap)) = 0, a system of 

three equations which cannot be solved for J (and hence to be discarded), or J = 0 , i.e. bi 2 = 621 . 
Hence, for a Kirchhoff-Love shell, ba /3 is necessarily symmetric (or in other words Aa /3 = A^a)- The 


three non-trivial strain incompatibility relations are 

*S'i 2 i 2 ~ ^ 11^22 + bi 2 = A'i 212 (^“, 0) and (4.2a) 

&al|2-&a2|l=A^12a3(0“,O). (4.2b) 

To rewrite the right hand side of the above relations in terms of dislocation density (torsion) we 

note that presently 

C) = 0 ) - C) gaAO^, C) 0 ) - C) 90.(0-, C) 0 ) (4.3) 

andC'3^(0“,O = T„^/3 r,O) (4.4) 

are the only non-zero components of the contortion tensor. Also, 

5 'a/3(6'", C) = aa/3(6'“) “ K + o(C) and gi3(0-, () = 5*3, (4.5) 

therefore, 

= a-^(e-) + 2Ca--a^%^(e-) + o(C) and g^^O-.Q = 5*3- (4.6) 

The strain incompatibility relations fl4.2l) can be rewritten as 

Si2i2 - 61,622 + 652 = -a„( 9 ‘>)|c 4 |, - C'12 + - C2",c;2]lj.„ (4.7a) 

and 6 (,i ,|2 - 6 (, 2 )|i = api( 6 “)C’J^| 3 ( 9 “. 0 ). (4.7b) 


With zero body force distribution and negligible inertia, the equilibrium equations for a Kirchhoff- 
Love shell take the form m- 


(cT'^“ + M^“ 6 ^)l„ + Mg“ 6 ^ = 0and (4.8a) 

(a^“ + M^X)bf,a - = 0, (4.8b) 
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where 


1 

2 


(4.9) 




I3a _ 


1 f dip 

2 V^ 


+ 


dpj \ 


and jM^'^ 


/ d^p ^^p \ 

dAgo^)' 


Equations fl4.7l) and fl4.8l) form the governing equations for the residual stress field and bending 
moment field for a Kirchhoff-Love shell with a continuous distribution of surface dislocation 
field specified by Ty^(6*",0). Dimensional analysis and representation theorem show that for a 
sufficiently thin isotropic Kirchhoff-Love shell, 'ip^Eais, Aa^) can be expressed as [22] 


= Eh(^C{J,, h) + h^Yl J2) 


i=3 


(4.10) 


where E is the Young’s modulus of the shell material and Ji := J 2 := , 

J3 := J4 := A^pA^,A-'^Af^\ J, := {E^pA^^A^^^A^^f, J, := A-\e^^ A^pE^,A^,A<^^ A^'^f 

and Jj := Ea^Ap^A^^AP'^A°'^^A^''. Here, C and Di are dimensionless functions. 


Remark 4-1 The components 0) of the torsion tensor do not contribute to the elastic defor¬ 

mation of a conventional Kirchhoff-Love shell. In other words, a Kirchhoff-Love shell geometrically 
admits only in-surface dislocation density represented by 0). This is also true for Kirchhoff- 

Love shells with uniform thickness distention, i.e. when d = [A + l)n with constant A. 


Remark 4-2 (Pure bending of an isotropic Kirchhoff-Love plate) In case of pure bending of a plate, 
= Aa 0 and = 0. The curvilinear coordinates ( 6 *^, 6 *^) can be identified with Cartesian co¬ 
ordinates; hence, = Aa 0 = Sa/s and Sa/S'y = 0 . Moreover, S '1212 = 0 . The strain incompatibility 
equations fl4.2l) are reduced to 


—A 11 A 22 + A 32 — A"i 212 (^°', 0 ) and 

(4.11a) 

~Ao- 1,2 + Ao- 2,1 = A'i2o-3(^'^, 0 ). 

(4.11b) 

The equilibrium equations fj4.8|) become [21] 


(d'^" + = 0 and 

(4.12a) 

(d^“ + = 0 , 

(4.12b) 

where cr“^ are to be interpreted as Lagrange multipliers cr"^(0") associated with the deformation 


constraint Ea /3 = 0 ; these are determined a posteriori after solving the complete boundary value 
problem. The bending moments are const it utively determined from the energy function 


'Ipi^afi)- 




1 P dip dip \ 

2 ^ 9Ajio/’ 


(4.13) 


where 


V’(Aa/3) = 


Eh^ 


24(1 


^AqqA^^ -j- (1 z/)Aq,^Aq,^ 


(4.14) 


u is the Poisson’s ratio of the shell material. 


18 
















Remark 4-3 (Small strain large rotation of an isotropic Kirchhoff-Love plate) Let and its 

spatial derivatives upto second order, be 0(e), where e := ^ << 1 and R is the minimnm principal 
radius of curvature of the shell mid-surface for a given deformation. Let Aq,^, and its spatial 
derivatives upto hrst order, be 0(e2). We identify (0^,0^) with the Cartesian coordinates on 
cu C Hence, = 5a0 and Hq/? = 0. A straightforward calculation shows that, upto 0(e), 

'Sq./3 • 2*^ (®cr/3,a A ^cra,P ~ A (^Eq-^ q^ E^jQ^ p Eq^p ^j^. (4.15) 

Consequently, the strain incompatibility equations for an isotropic Kirchhoff-Love plate take the 
form (upto 0(e)) 


2 £'i 2 ,i 2 ~ -E' 11,22 ~ -E' 22,11 ~ A 11 A 22 + A ^2 — A'i 212 (^“, 0) and (4.16a) 

— Ao-1,2 + Ao-2,1 = A'i2o-3 (6*“, 0 ). (4.16b) 

Moreover, Ji = 0(e), J 2 = O(e^), J 3 = 0(e), J 4 = 0(e), J 5 = 0(e^'^^), Jq = 0(e^'^^), and 

J 7 = O(e^). The quadratic strain energy function, upto O(e^), neglecting the non-conventional 
coupling term J-j (it is important to note that, unlike the small deformation theory where the 
strain energy is decoupled at O(e^), a coupling term is always present), is 


'0(T/a^, Aq;^) 


Eh 


2(1 - z/2) 


^Efy^Q^Epp (1 Eo^pEoip I -|- 


Eh^ 


24(1 - 1^2 


^E(^aA.pp (1 h'')A.apEap 

(4.17) 

The corresponding equilibrium equation in terms of E^p and Kap can be written upto 0(e) using 
fl4.8p . fl4.15l) . and the above energy function. 


Acknowledgement: We express our gratitude to Prof. Robert Bryant for the insightful discussion 
at mathoverflow.net [2] about existence and uniqueness issues of the non-linear PDEs encountered 
in the proof of Theorem 3.2. 
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